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The Gibbs-free-energy landscape for the solute association in nanoconfined aqueous solutions* 
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The theoretical model and the numerical analyses on the Gibbs-free-energy of the association states of am- 
phiphilic molecules in nanoconfined aqueous solutions are presented in detail. We exhibit the continuous change 
of the Gibbs-free-energy trend, which plays a critical role in the association states of the system transforming 
from the dispersion state, through the “reversible state”, and finally to the aggregation state in amphiphilic 
molecule solutions. Furthermore, for the “reversible state”, we present the difference in the free-energy bar- 
rier heights of the dispersion state and aggregation state, resulting from the competition between the entropy, 
which makes the solute molecules evenly disperse in the solution and the energy contribution driving the am- 
phiphilic molecules to aggregate into a larger cluster. These findings provide a comprehensive understanding 
of confinement effects on the solute association processes in aqueous solutions and may further improve the 


techniques of material fabrication. 
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I. INTRODUCTION 


The solute association in aqueous solutions plays a funda- 
mental role in various processes, including crystal growth [1, 
2], colloid particles nucleation [3, 4], surfactant micelliza- 
tion [5, 6], gas hydrate formation [7, 8], and biomolecule re- 
lated applications [9-11]. When the confined environment, 
such as the space between the macromolecules in a cell [12], 
pores in the soil [13], and nanochannels [14, 15] in man- 
made nanomaterial, is involved, the solute molecules exhibit 
different behavior from that in bulk. Zhang et al. [16] found 
that the stability of the native state of the protein can be en- 
hanced at a smaller confinement volume. For the surfac- 
tant solutions confined in the nanochannel, various morpholo- 
gies during the transition appear [17]. The experiment in the 
confined geometry designed by Julie et al. [18] reported that 
a single skeleton crystal is observed for sodium chloride at a 
high supersaturation rate. 

Our previous work [19] has shown that the pentanol 
molecules can switch between the dispersion state and the 
aggregation state, called “reversible state transition’, in 
nanoconfined aqueous solutions with the aid of molecular dy- 
namics simulations. We built the ideal physical model and 
numerically calculated the Gibbs-free-energy of a solute with 
respect to the number of molecules in the cluster. It is found 
that the unexpected observation is attributed to the two min- 
ima separated by a maximum in the Gibbs-free-energy curve 
under the confinement and existence of a free-energy bar- 
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rier, which is the same order as the thermal fluctuations. The 
confinement also leads to the enhancement of the critical ag- 
gregation concentration, which results from the fact that the 
solute number is not sufficient to sustain the stable cluster at 
a lower concentration in the confined volume. 

The pentanol molecule in our MD simulation is described 
by the Gromos force field, where either a methyl group or 
ethyl group is considered as a whole group and their hydrogen 
atoms are implicitly included. This way, the acceleration of 
computation is evident to track the dynamic association pro- 
cesses, including the dispersion, reversible switch, and aggre- 
gation states. The packing of pentanol into clusters is of little 
influence, as these groups always behave like a whole entity 
during the aggregation processes, although a rough descrip- 
tion of water structures around the methy] or ethyl groups will 
be obtained because hydrogen atoms on these groups are not 
explicitly considered. 

In this study, we first derive the free-energy formula in 
terms of the ideal physical model in detail to illustrate the 
landscape of Gibbs-free-energy of amphiphile association be- 
haviors in the nanoconfined geometry. Secondly, especially 
for the “reversible state”, we discuss Gibbs-free-energy barri- 
ers for the dispersion and aggregation states and numerically 
investigated the relationship between the Gibbs-free-energy 
barrier height and the concentration of solute molecules. Fi- 
nally, the conclusion is summarized based on the characteris- 
tic of the “reversible state” of the amphiphilic solute associa- 
tion behavior in nanoconfined aqueous solutions. 


Il. MODEL AND METHODS 


A. Model of the amphiphilic molecule and cluster 


(1) A single amphiphilic alcohol molecule is modeled by a 
cuboid with a hydrophilic head and a hydrophobic tail, 
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Fig. 1. (Color online) (a) Molecular models for the single am- 
phiphilic molecule and (b) the cluster with a radius R. The cuboids 
in red and green represent the hydrophilic head (-CH,-OH) and hy- 
drophobic tail (alkyl chain). The molecule width and length and the 
hydrophobic tail length are denoted by a, L, and l, respectively. In 
the cluster model, the densely packed hydrophobic tails, shown in 
green, prefer to stay in the interior core and the hydrophilic heads, 
shown in red, are on the surface. The surrounding water environment 
is not shown for clarity. 


as illustrated by the red and green geometry in Fig. 1(a). 
The molecule width, molecule length, and tail length are 
denoted by a, L and l, respectively. The volume of a 
molecule is Ym = a? L. 
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where kg is the Boltzmann constant and T is the system tem- 
perature. The first three terms in the curly brace are con- 
tributions from the surface energy of hydrophobic tails and 
the fourth term is from placing the hydrophilic heads onto 
the cluster surface. The n - |Ajitransfer| is the free energy 
difference of transferring a single hydrophobic tail from the 
dispersed phase into the hydrophobic parts of the cluster. The 
terms in the square brackets represent the free-energy contri- 
bution from the N — n dispersed molecules. Ypyoo and Yyoo 
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(2) A cluster consisting of densely packed amphiphilic 
molecules in water is treated as an ideal spherical aggre- 
gate [20] with the radius R(n): 


R(n) = (3num/4r)! = (3na?L/4r)*/3, (1) 


where n is the number of molecules in the cluster, namely, 
the aggregation number of the cluster. The hydrophobic 
tails prefer to stay in the interior core and the hydrophilic 
heads like to occupy the interfacial region exposed to the 
surrounding water environment, as shown in Fig. 1(b). 


B. Methods for Gibbs-free-energy calculations 


The theoretical description of Gibbs-free-energy is in- 
troduced to calculate the free-energy of amphiphilic solute 
molecules in nanoconfined aqueous solutions. We consid- 
ered that only one single cluster is formed and the solute 
molecule is either in the cluster or dispersed in the solution. 
The total amount of solute is denoted by N and the number 
of molecules in dispersion is N — n. The Gibbs-free-energy, 
G, as the function of both the cluster aggregation number, n, 
and the solute number JN, is defined as 


al 
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are the interfacial tension coefficients of hydrophobic tails- 
water and water-vapor and ôhw, dy are the corresponding Tol- 
man lengths [20]. The confinement effect is displayed by V, 
which is the volume, excluding the walls. 


The relationship between the aggregation number, n, of the 
cluster and the amount of solute, N, namely, n( N), can be 
derived by differentiating the Gibbs-free-energy, G (n, N), to 
n. The n(N) is given by the relation 


1 
4 3Um \ 3 
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where the terms on the both sides imply the chemical bal- 
ance between the cluster and other dispersed molecules in the 
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confined aqueous solutions. By defining the solute concentra- 
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tion, C, as C = N/V, other expressions, G(n, C) and n(C), 
can be obtained from G(n, N) and n(N) at a certain volume. 

To proceed with the numerical calculations, the 
coefficients are parameterized as follows: T = 300K; 
V = 5.0nm x 5.0nm x 5.0nm; For the water, yoo = 72 mN/m 
and 6, = 1.23 Å; For the pentanol molecule, Ynvoo = 49.53 mN/m, 
a = 2.60A,1 = 5.97A, b = 0.02, and —|Aptuansfer] — B = 
—5.26ksT. All of these numerical analyses are achieved by using 
the Maple 18 software. 

With respect to the aqueous solutions solvated with the amount 
of amphiphilic solute, N, in a certain volume, V, the Gibbs-free- 
energy, G(n, N) or G(n, C), can be regarded as the function of the 
aggregation number of the cluster n and thus it is denoted by G(n) 
for clarity. The local minimum in the Gibbs-free-energy, G(n), is 
found by 


dG(n) d?G(n) 
= d : 4 
aa O an In? <0 (4) 
Similarly, the local maximum is found by 
dG(n) d?G(n) 
= d ; 
dn QO an In? >0 (5) 


The barrier height, Garier, is defined by the magnitude of the 
difference of Gibbs-free-energy values between the local minimum 
and its adjacent maximum, i.e., 


Gharrier = |G max = Gwin. (6) 
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Fig. 2. (Color online) The mathematical picture for the relation, 
n(C), formulating the dependence of the aggregation number on 
the cluster, n, on the solute concentration, C. 


Il. RESULTS 


To describe how many solute molecules will aggregate into a clus- 
ter and form an association state of the system as the solute concen- 
tration increases, we plot the n(C) function in terms of Eq. (3). As 
shown in Fig. 2, three regions of solute concentration, C’, can be rec- 
ognized. When C < 0.266 mol/L, the aggregation number is n ~ 1, 
which indicates no aggregates formed and the solute molecules like 
to disperse in the solution. When sufficient solute molecules are 
added into the solution, i.e., C > 0.425 mol/L, an evident cluster ap- 
pears and it grows larger as the solute concentration increases. Inter- 
estingly, when 0.266 mol/L < C < 0.425 mol/L in the “reversible 
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Fig. 3. (Color online) (a) The Gibbs-free-energy, G (n), of the so- 
lute molecules in confined aqueous solutions versus the aggregation 
number of the cluster, n, under various solute concentrations. The 
subpicture shows the G(n) when C = 0.279 mol/L. The G(n) 
curves are shifted for easy comparison. (b) The landscape of the 
Gibbs-free-energy with respect to the aggregation number n and so- 
lute concentration, C’. The surface is colored according to the value 
of the Gibbs-free-energy. 


state”, the solute can exhibit both the dispersion state and the aggre- 
gation state, because the three aggregation numbers, n, correspond 
to a certain solute concentration. This can be seen more clearly in 
the Gibbs-free-energy landscape in Fig. 3. 

In terms of the Gibb-free-energy formula in Eq. (2) and gradu- 
ally increasing the solute concentration C, the images of G(n, C) 
at various C are shown in Fig. 3. For C = 0.266 mol/L, there is 
only one minimum where the aggregation number, n ~ 1, corre- 
sponding to the dispersion state. For the system with a larger solute 
concentration of C = 0.439 mol/L, the minimum at n ~ 26 demon- 
strates that most of the molecules aggregate together, thus forming 
a cluster. However, when C = 0.279 mol/L, the G(n) displays two 
minima at n ~ 1, 12 separated by a maximum at n ~ 4. According 
to the definitions given by Eq. (4) and Eq. (5), these values corre- 
spond to the stable dispersion state, the stable aggregation state, and 
the metastable state, respectively. This bistability also appears when 
C = 0.332 mol/L and 0.385 mol/L. The continuous change of the 
Gibbs-free-energy trend with increased solute concentration is more 
clearly shown in the landscape of the Gibbs-free-energy. Figure 3(b) 
presents the Gibbs-free-energy landscape with respect to the aggre- 
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Fig. 4. (Color online) The dependence of the dispersion barrier 
height and the aggregation barrier height on the solute concentra- 
tion, C, when the system is in the reversible state. 


gation number and the solute concentration. Overall, as the solute 
concentration, C, increases, the G'(n) exhibits the continuous vari- 
ation that the minimum corresponding to the dispersion state disap- 
pears when accompanied by the appearance of the aggregation state. 

The emergence of the minima and maximum for the solute con- 
centration C = 0.279, 0.332 and 0.385 mol/L in Fig. 3 makes it 
possible for the system to switch between the dispersion state and 
the aggregation state by overcoming the free-energy barrier under 
thermal fluctuations. This activation process at a certain temperature 
directly relies on the height of barrier. Considering the bistability in 
the “reversible state” region 0.266 mol/L < C < 0.425 mol/L, the 
barrier heights for the dispersion state and the aggregation state ver- 
sus the solute concentration, C’, are presented in Fig. 4 by following 
the definition of the free-energy barrier, Garnier, in Eq. (6). 

The max barrier height is the order of hydrogen bond energy 
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(~ 8ksT) in water at 300K. That is to say, the transition of so- 
lute molecules can occur between these two different states under 
the thermal fluctuations. Because the energy contribution driving 
the aggregation dominates the behavior of the solute and the larger 
cluster owns the lower free-energy, it is seen that the barrier height 
for the aggregation state rises with the increase of the solute con- 
centration, C. The disappearance of the dispersion energy barrier is 
probably attributed to the fact that the entropy effect promoting the 
even distribution of the solute in the solution is reduced. 


IV. CONCLUSION 


In summary, we presented the detailed derivation of the Gibbs- 
free-energy landscape and describe the amphiphilic solute associ- 
ation behavior in nanoconfined aqueous solutions. Numerical cal- 
culations showed that under the confined geometry, as the solute 
concentration increases, the solute free-energy curve shows the con- 
tinuous transition from one minimum with a aggregation number of 
n ~ 1, through two minima separated by one maximum, finally to 
one minimum where n ~ N. This change indicates the system will 
transit from the stable dispersion state, through the “reversible state”, 
and finally to the stable aggregation state. In the “reversible state” 
region, the max energy barrier is the order of the hydrogen bond 
energy (~ 8ksT) in water, which allows the system to switch be- 
tween the dispersion state and the aggregation state. The difference 
between the free-energy barrier for the dispersion state and the ag- 
gregation state is identified. The decrease in the dispersion barrier 
height and the rise of the aggregation barrier height resulted from 
the competition between the entropy making the solute evenly dis- 
tributed in the solutions and the energy contribution driving the so- 
lute aggregate into the cluster, as the solute concentration increases. 

Our work provides a comprehensive understanding of the 
confinement on the behavior of solute molecules in aqueous solu- 
tions. It is expected that the theory can also be applied to other 
amphiphilic molecules with more complicated structures, such as 
peptides or surfactants, in a nanoporous medium or cell, which are 
widely spread confined environments. 
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